clear all; clc; close all
alfa_f = deg2rad(10);
beta_f = deg2rad(0);

alfa_s = deg2rad(-90:5:90);
beta_s = deg2rad(-90:5:90);
beta = H(alfa_f,beta_f);
for a=1:length(alfa_s)
    for b=1:length(beta_s)
        S = 1;
        y=H(alfa_s(a),beta_s(b))*S;
        U = beta' * y;
        F(a,b)= abs(U/S)^2;
        
    end
end


[alfa_sm,beta_sm] = meshgrid(alfa_s,beta_s);
figure;
surf(alfa_s,beta_s,F)
figure;
[Fx,Fy,Fz] = sph2cart(beta_sm, alfa_sm, F);
surf(Fx,Fy,Fz)
